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Complex numbers are an intrinsic part of the mathematical formalism of quantum theory, and 
are perhaps its most characteristic feature. In this paper, we show that the complex nature of 

the quantum formalism can be derived directly from the assumption that a pair of real numbers is 
associated with each sequence of measurement outcomes, with the probability of this sequence being 
a real- valued function of this number pair. By making use of elementary symmetry conditions, and 
without assuming that these real number pairs have any other algebraic structure, we show that 
these pairs must be manipulated according to the rules of complex arithmetic. We demonstrate 
that these complex numbers combine according to Feynman's sum and product rules, with the 
modulus-squared yielding the probability of a sequence of outcomes. 



I. INTRODUCTION 

Complex numbers are perhaps the most characteristic mathematical feature of quantum theory. In recent years, 
there has been growing interest in elucidating the physical origin of this (and other) mathematical features of quantum 
theory by deriving - or reconstructing - the quantum formalism from one or more physical principles, and signicant 
progress has been made [1-19]. For example, many approaches are able to derive specific equations or predictions 
such as Schroedinger's equation or Malus' law [5, 6, 8, 9, 17, 18]. However, the derivation of a significant part of the 
quantum formalism has thus far relied upon abstract assumptions such as the introduction of the complex number 
field [10, 12] or upon several disparate features of quantum phenomena [13 15]. 

In this paper, we present a novel reconstruction of Feynman's reformulation of quantum theory [20] . Our approach 
differs from previous approaches in that it avoids ad hoc introduction of the complex number field, and in that it rests 
essentially on a single postulate, namely: 

Pair Postulate: each sequence of measurement outcomes obtained in a given experiment is represented 
by a pair of real numbers, where the probability associated with this sequence is a continuous, non-trivial 

function of both components of this real number pair. 

Since the probability is the only information which is accessible in a given experiment, this postulate expresses the 
simple idea that it requires twice as many degrees of freedom to describe a physical system than one can access through 
a given measurement. This idea has played an important role in some previous attempts to reconstruct quantum 
theory [40], and can also be regarded as one way of stating Bohr's principle of complementarity [21]. 

Using symmetry and consistency conditions that arise naturally in an operational framework, and making a few 
elementary physical assumptions, we show that this postulate leads to Feynman's riilcs of quantum theory. Most 
importantly, we show that the number pairs assigned to each sequence of measurement outcomes must be manipulated 
according to the rules of complex arithmetic, without assuming this at the outset. Specifically, in the language in 
which Feynman's rules are usually expressed [20], we show that, if the pair associated with a path which a system 
classically can take from an initial event, Ei, to a final event, Ef, is written as a complex number, or amplitude, then: 

(a) If a system classically can take more than one path from Ei to Ef, then the total amplitude for the transition 
is given by the sum of the amplitudes associated with these paths. 
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(b) If the transition from Ei to Ef takes place via intermediate event Em, the total amplitude is given by the product 
of the amplitudes for the transitions Ei — > Em and Em —>£"/, and 

(c) The probability of the transition from E^ to Ef is proportional to the modulus- squared of the total amplitude 
for the transition. 

Our approach is partly inspired by two previous reconstructions of Feynman's rules due to Tikochinsky [7] and 
Caticha [11]. Tikochinsky postulates that a complex number is associated with each path that a system can take 
between events. By identifying a set of symmetries associated with these paths, he adapts an argument used by Cox 
to derive probability theory [22, 23] to show that these complex numbers must combine according to Feynman's rules. 
Caticha's approach is similar, except that he operationalizes the classical notion of 'path' within an experimental 
framework. Both of these authors assume at the outset that complex arithmetic is to be used, and also implicitly 
assume that certain complex functions are analytic. Such assumptions are given no a priori physical justification, 
which detracts from the physical insight that these reconstructions can provide. We show the complex structure of 
quantum theory need not be assumed, but can in fact be derived. The remainder of this paper is organized as follows. 

In Section II, we present an experimental framework which provides the basis for the reconstruction. The exper- 
imental framework provides a fully operational language which we use in place of the classical language of 'paths' 
employed in Feynman's original formulation. Sequences of measurement outcomes, and two operators that can be 
used to combine them in series and in parallel, are introduced. Five fundamental symmetries associated with these 
operators are derived. We then obtain a representation of this space of sequences by representing each sequence by a 
number pair and by requiring that these pairs combine through pair operators which share the same symmetries. 

In Section III, we use the symmetry constraints on the pair operators to determine their form. This restricts 
their form to a few possibilities. We then impose the requirement that the probability associated with a sequence 
is determined by the number pair that represents that sequence. This requirement eliminates all but one of these 
possibilities, and completes the derivation by yielding the modulus-squared relation between the probability and the 
pair. 

We conclude in Section IV with a discussion of the results obtained and of potential future developments. 

II. EXPERIMENTAL FRAMEWORK 

We consider experimental set-ups in which a physical system is subject to successive measurements Mi, M2, M3, . . . 
at successive times fi, f2, fs, . . . . The system is allowed to undergo interactions in the intervening intervals. We 
summarize the outcomes obtained in a given run of the experiment as a sequence A = [mi, m2, ma, . . . ]. The 
measurements can be of different features of the system, but we shall label the outcomes of each measurement 
as 1, 2, 3, 4, . . . as far as needed in each case. 

Consider, for example, the Stern-Gerlach set-up shown in Fig. 1. Here, a source supplies silver atoms which pass 
through the apparatus, undergoing successive measurements of components of spin. Each measurement is performed 
by a magnet equipped with two wire-loop detectors (as sketched in the figure) which do not absorb the atoms. Between 
the measurements, the spins may interact with a uniform magnetic field. For silver atoms, it is found experimentally 
that each measurement can only have two possible outcomes, which we label 1 and 2. These measurements are 
repeatable in that the same result is always obtained if the same measurement is immediately repeated. 

We might, for example, obtain the sequence A = [2,1,2], or perhaps B = [2,2,1]. Under repeated trials of this 
experiment, the probability distribution over the outcome of M3 is observed to be independent of any interactions 
the system had prior to M2, including the outcome of Mi. In such a case, we say that the earlier measurement M2 
establishes closure with respect to the later M3 [41]. Closure, in which current information overrides past information, 
is a basic feature of experiments on quantum systems. 

We can also set up coarser experiments, such as the one shown in Fig. 2. Here, the measurement M2 performed at t2 
uses only a single detector whose field of sensitivity includes outcomes 1 and 2 of M2 in the original experiment. Now, 
if the coarser M2 registers an atom, only outcome 1 or 2 could be obtained if measurement M2 was then performed 
immediately afterwards. Accordingly, we write the outcome of M2 as (1,2), and we say that the measurement M2 
coarsens outcomes 1 and 2 of the original M2. Using M2, outcome (1,2) can be refined to finer outcomes 1 and 2, 
but those latter outcomes cannot be further refined. An outcome that cannot be further refined is said to be atomic. 
If a measurement, such as M2, has all of its outcomes atomic, we shall call the measurement itself atomic. The 
notation for non-atomic outcomes is naturally extended to the case where an outcome can be refined into more than 
two outcomes. 

Generalizing the Stern-Gerlach example, we consider set-ups where the measurements (of a particular property) are 
repeatable and either atomic or coarsened versions of such, and where the first and last measurements in the set-ups 
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Ml Ma M3 

FIG. 1: Schematic representation of a Stern-Gerlach experiment performed on silver atoms. A silver atom from a source (an 
evaporator) is subject to a sequence of measurements, each of which yields one of two possible outcomes registered by non- 
absorbing wire-loop detectors. A run of the experiment yields outcomes mi , m2 , ms of the measurements Mi , M2 , M3 performed 
at times ti,t2,t3, respectively. 

are atomic. We also take the observed system to be sufficiently simple that the atomic measurements establishes 
closure with respect to any future measurement, and that any interaction with the system between measurements 
preserves this closure. 

A. Combining Sequences 

We now consider different ways in which sequences of measurement outcomes can be combined with one another to 
generate other sequences. We use two kinds of relations between sequences, namely parallel and series combination. 

1. Sequences in Parallel 

First, consider an experimental set-up consisting of three measurements. Mi, M2 and M3 performed in succession. 
On one run, this generates sequence A = [mi, 1712,111,3] and, on another run, sequence B = [mi,m'2,m3], with m2 ^ 
m'2. Then consider a second set-up, identical to the first except that the intermediate measurement M2 coarsens 
outcomes m2 and m'2 of M2, and suppose that this generates the sequence C = [mi, (7712,7722), ms]. We shall say that 




Ml M2 M3 



FIG. 2: A Stern-Gerlach experiment where the field of sensitivity of the intermediate measurement detector spans the fields of 
sensitivity of both of the detectors of the corresponding measurement in Fig. 1. 
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the sequence C combines A and B in parallel (Fig. 3). We symbolize this relation by defining a binary operator, V, 
which here acts on A and B to generate the sequence 

C = AyB. (1) 

Generally, the binary operator V combines any two sequences obtained from the same experimental set-up differing 
in only one outcome. 




A y B = C 

FIG. 3: Combination of sequences in parallel. Graphical depiction of the sequences A = [mi,m2,m3], B = [mi, m2, ma], 
and C = [mi, (m2,m2),m3], respectively. 

From the above definition, it follows at once that V is commutative and associative. To establish the first, notice 
that 

BV A= [mi, {1712,1712), m^], 

and since (rn2,m2) = (m2,TO2), it follows that V is commutative, 

AV B = BV A. (2) 

To establish the second property, consider the three sequences A = [mi,m2,ms], B = [mi,m'2,ms\, and C = 
[mi, 70.2,1713], with 1712, m'2 and m2 distinct. 
These sequences can be combined to form D = [mi, (m2, mg, 1112), m^] in two different ways, namely 

D = (AVB)VC and L' = AV(BVC), 

which implies that V is associative, 

(^VB)VC = ^V(BVC). (3) 

2. Sequences in Series 

Consider the two sequences A — [mi, m2] and B — [m2, m^], in which outcome m2 is the same in each (see Fig. 4). 
We now define the binary operator, • , which chains two such sequences in series. This acts on A and B to generate 
the sequence 

C = A- B = [mi, 702,70,3]. (4) 

Generally, the binary operator • combines together any two sequences obtained from experimental set-ups where the 
last measurement (and the outcome) of one sequence coincides with the first measurement (and the outcome) of the 
other. 

By considering the three sequences A = [mi, m2], B = [m2, m^], and C = [ma, m4], we see that • is associative, 

{A-B)-C = A-{B-C). (5) 
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A • B = C 

FIG. 4: Combination of sequences in series. Graphical depiction of sequences = [mi,m2],S = [m2,m3] and C = [mi,m2,m3], 
respectively. 

Finally, consider the sequences A = [1711,1712, TO3] and B = [mi, mj, TO3] and C = [ma, 7174]. These can be combined in 
two equivalent ways to yield D = [mi, (m2, 1712), m^, 1714], namely 

D = {AVB)-C and D = {A-C)V (B-C). 

Hence, the operation • is right-distributive over V, 

{Ay B)-C = {A-C)y {B-C). (6) 

Similar considerations show that • is also left-distributive over V, 

C-{AyB) = {C-A)y{C-B). (7) 



B. Sequence Pairs 

Following the first part of our pair postulate, we represent each sequence, A, by a real number pair, a = (ai, 02)^. 
We have determined that the parallel and series operators, V and •, possess the symmetries given in Eqs. (2), (3), 
(5), (6), and (7). These symmetries must be reflected in the representation. For example, if pairs a, b represent the 
sequences A, B, respectively, then the pair c that represents C = Ay B must be determined by a, b through the 
relation 



c = ae b. 



(8) 



where ® is a pair-valued binary operator, assumed continuous, to be determined. Then, since V is commutative, we 
have Ay B = By A, so that 

aeb = bea. (SI) 

In addition, since V is associative, 

(a®b)ec = a® (b®c). (S2) 

Similarly, if the sequences A,B and C are related by C = A-B, then the pair c that represents C must be 
determined by a, b through the relation 

c = a b, (9) 

where © is another pair-valued binary operator, assumed continuous, also to be determined. Prom the associativity 
of • , it follows that also has associative symmetry, 

(a0b)0c = a0(b0c). (S3) 

Finally, since • is right- and left-distributive over V, it follows that the pair operators also have distributive symmetry, 

(a®b)0c= (a0c)® (b0c) (S4) 

a0 (b® c) = (a® b) ® (a0c). (S5) 
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III. DERIVATION OF FEYNMAN'S RULES 

In Sec. Ill A, we shall use the symmetry equations (S1)-(S5) to fix the form of ® and to restrict to one of 

five possible forms. Then, in Sec. IIIB, wc shall introduce a connection between pairs and probabilities. This will 
restrict to a unique form, and fix the functional connection between pairs and probabilities. 

A. Solution of the Symmetry Equations for ® and 

1. Solution of Commutativity and Associativity Equations for ® 

Commutativity and associativity of impose strong constraints on the possible forms that the operator can take. 

To illustrate the nature of this constraint, consider a binary operator o which acts over the real numbers. In this 
one-dimensional case, there exist a number of theorems which show that, if operator o is continuous and associative, 
and possesses a small number of additional properties [42] , then the operator must satisfy the equation 

f{xoy) = f{x) + f{y), (10) 

where / is a continuous and strictly monotonic function. That is, given a binary operator over the reals satisfying the 
above-mentioned conditions, one can always invertibly transform the real-line such that, in the transformed space, the 
operator o is represented by the addition operator. Hence, without any loss of generality, one can choose to perform 
the composition operation in the transformed space. Parenthetically, this result forms the basis of Cox's derivation 
of probability theory [22, 23] and is the rationale for additivity in measure theory [24]. 

In the two-dimensional case with which we are concerned here, an analogous result holds, namely that, for contin- 
uous, associative and commutative 0, then the operator must satisfy the equation 

F(a©b) = F(a)-hF(b), (11) 

where F is an invertible and continuous pair- valued function. This was proved by Aczel and Hosszii [25] with the aid 
of minor technical assumptions. 

Hence, without any loss of generality, we can transform the space of pairs such that the operator becomes 
represented in standard form by the additive operator. Explicitly, in this standard form, 

which we refer to as the sum rule. Note that the only freedom left in the sum rule is a real invertible linear 
transformation of the space of pairs, 

ill) = {u v) ill) ' 

with SV — TU ^ 0. That is, whenever the sum rule holds between pairs, it also holds between the corresponding 
transformed pairs. We shall make use of this fact below. 

2. Solution of Associativity and Distributivity Equations for Q 

Having shown that corresponds to component-wise addition of number pairs, we proceed to show that corre- 
sponds to a form of multiplication. 
a. Distributivity of (D . First, define 

G(a,b) = a0b, 

where the pair- valued function G is to be determined through Eqs. (S4) and (S5), which become 

G(a -h b, c) = G(a, c) + G(b, c) 
G(a, b c) = G(a, b) + G(a, c). 



(13) 
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Defining ra = r(ai,a2)^ = {rai,ra2)^ , in accordance with Eq. (12), with r real, it follows that 

G(ria,r2b) = rir2G(a,b). 

Introducing two-dimensional basis pairs ei and 62, it then follows that 

G(a, b) = G(aiei + 0262, foiei + 6262) 

= ai6iG(ei,ei) + oi62G(ei,e2) + 026iG(e2,ei) + a262G(e2,e2) 

ft) 

where 7 = (71, 72, 73, 74; 75, 76, 77, 78) is a real-valued vector to be determined, in which the semicolon partitions 7 into 
components that, respectively, effect the first and second part of the real pair. Hence, the left- and right-distributivity 
of over ® implies that a b has the bilinear multiplicative form 

('^^^ & ( = ( ^'^^^^'^ T2"l^2 + 730261 + 74a2&2\ 

\a2j ' \b2j ^7501^1 +7601^2 +77a2&i +7801262/ ■ ^ ^ 

b. Associativity of Q. Substituting this form of a0b into the 0-associativity condition, Eq. (S3), and solving the 
resulting equations (see Appendix A), one finds that 7 can take one of three possible forms, namely a commutative 
form 

7 = (6» - Ve, 0e, 0e, (j); 0e, 0, 0, ip + (j)e), (15) 

with real constants 6, (j), tp, e, and two non-commutative forms 

7 = {0, <f>, 0, 0; 0, 0, e, <P) (16) 
7 = {0, 0, ^, 0; 0, 0, 0, V). (17) 

Using the freedom described in Eq. (13), we can transform these solutions to standard forms. To do so, we note 
that, under the transformation of Eq. (13), the relation c = a0 b transforms to 

/ c[\ ^ fS T\ f 7iai6i + 720162 + 73a26i + 74a262A 
\P V V^sfli^i +750162 + 770261 +730262/ 

^ ( 7Wi6i + 72ai62 + 73«26'i + iia'^b'^\ 
yi^a'^h'^ + 76«i6'2 + 77a26'i + 78462^ ' 

where 

and where 7' = (7^, ... ,73) is the representation of 7 in the space of the transformed pairs. Equating coefficients 
of a'l, 6'i, 62 identifies 7' as 
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(19) 



Using this transformation, Eqs. (15), (16) and (17) can be reduced to standard forms. 
In particular, Eq. (15) takes the standard form 



7 =(1,0, 0, 11; 0, 1, 1,0), 



(20) 
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where /i = sgn{4d(j) + ip"^) can be —1, 0, or +1. Through Eq. (14), case ^ = —1 gives what we recognize as complex 
multiplication, while the cases /z = and /x = +1 give variations thereof. The transformation needed to recover 
Eq. (15) from this standard form is 

S T\ _1 [20- V'e 2(l)e + tp' 
U V) ~ 2y eA 

where 

^^fvW+V^ ifM = ±l 
[1 if/x = 0. 

When SV — TU ^ 0, the inverse of this transformation exists, so that Eq. (15) can be returned to the standard form, 
Eq. (20). Note that this standard form with /x = +1 can be reached from the even simpler form 

7=(1, 0, 0, 0;0, 0, 0, 1), (21) 

that we use later, by applying the invertible transformation 

S T\ _(\ -1 
V V)~\\ 1 

When, on the other hand, SV — TU — 0, the transformation would be singular, hence disallowed. This would 
happen if ^ = where ( = ip + (j)€. Eq. (15) would then have been 

7 = (0e^, (^^e, (pe, (p; Ce^, Ce, (e, C), 

Observing the linear relation ^Ci = (pC2 between the components of any product c = a© b thus defined, we note that 
this could be transformed by rotation to c' = ((Z^, 0)^. This lacks the two components that we demand of an arbitrary 
pair, so the singular case is inadmissible. 

Continuing in this style, Eq. (16) takes the standard form 

7 =(1,0, 0,0; 0,0, 1,0). (22) 

The transformation needed to recover Eq. (16) from this standard form is 

S T\ _ f (P 
U VJ ~ \-(p 

which is invertible unless and (p both vanish. Similarly, the other non-commutative form Eq. (17) has the standard 
form 



This transforms to Eq. (17) through 



7 =(1,0, 0,0; 0,1, 0,0). (23) 



S T\ _ f iP 
U V) ~ \-ip 



which is, again, invertible unless and ip both vanish. 

In summary, imposing associativity of restricts 7 to one of five possible standard forms, 

7 = (1, 0, 0, -1;0, 1, 1, 0) (CI) 
7 = (1, 0, 0, 0; 0, 1, 1, 0) (C2) 
7 = (1, 0, 0, 0; 0, 0, 0, 1), (C3) 



and 



7 = (1, 0, 0, 0; 0, 1, 0, 0) (Nl) 
7 = (1, 0, 0, 0; 0, 0, 1, 0). (N2) 



each of which, through Eq. (14), defines a way to multiply pairs. The first three give complex multiplication (CI) 
followed by two variations thereof (C2 and C3), and the last two give non-commutative multiplication (Nl and N2). 
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B. Probability of a Sequence 

At this point, symmetry alone can take us no further in determining the precise form of the operator 0. In order 
to make progress, wc make use of the second part of our pair postulate, and introduce a connection between the pair 
that represents a sequence and the probability associated with that same sequence. 

We define the probability P{A) associated with sequence A = [mi,m2, . . . ,m„] as the probability of obtaining 
outcomes m2, . . . , m„ conditional upon obtaining mi, 

=Pr(m„,m„_i,...,m2|mi). (24) 

Following our pair postulate, we now require that P{A) is determined by the pair, a, that represents sequence A, so 
that, for any a, 

P(A)=p(a), (25) 

where p is a continuous real- valued function that depends non-trivially on both real components of its argument [43] . 

Our goal in this section is to determine the constraints imposed by probability theory on the form of p and, in the 
process of doing so, to show that only form (CI) can yield a form of p which meets our stated requirements. 



1. Probability Equation 

Consider the two sequences A = [mi,m2] and B = [m^, m^] of atomic outcomes. Since outcome m2 is the same in 
each, C = A - B ]s given by C = [m\,m2,mz\- The probability, -P(C), associated with sequence C is given by 

P{C) = Pr(m3, 7712 I mi), 

which, by the product rule of probability theory, can be rewritten as 

P{C) = Pr(m3 I m2, mi) Pr(m2 | mi). 

Since m2 is atomic, measurement M2 (with outcome m2) establishes closure with respect to M3 (with outcome ms), 
by definition overriding the earlier outcome mi. Therefore, the probability of outcome ms is independent of mi, and 

the above equation simpUfies to 

P{C) = Pr(m3 I m2) Pr(m2 | mi) 
= P{B)P{A). 

Hence, for any a, b, the function p must satisfy the equation 

p(a0b) =p(a)p(b). (26) 
Solving for the function p that satisfies this equation in each of the five forms of 7 given above, we obtain 

Case CI: p(a) = {al + alY'^; 

Case C2: p(a) = |oi|"e''"2/'^i; 

Case C3: p(a) = \ai\'^\a2\^ ; 

Case Nl: 39(a) = |ai|"; 

Case N2: p{a) ^ |ai|"; 

with Q:, d real constants (see Appendix B). The solutions for p in the case of the two non-commutative forms (Nl) 
and (N2) depend only on the first component of its argument. That is not admissible, so those two forms are rejected. 
Of the five possible forms of 7, we are left with three: (CI), (C2) and (C3). 
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2. Reciprocity 

Suppose that the sequence A = [mi,n2\ is obtained from an experiment where measurements M and N are 
performed at times ti and t2, respectively. Now consider the experiment where the measurements are performed in 
the reverse order, so that N is performed at time ti, followed by M at time ^2, and suppose that the sequence obtained 
is A = [n2,mi], where the over-arrow symbolizes a unary operator acting on the sequence A. 

Suppose M and N arc Stcrn-Gerlach measurements as in Sec. II. Then, in the limit as t2 ^ ti, it follows from 
rotational symmetry and the empirical fact that the outcome probabilities of N depend only on the magnitude of 
the angle between the axes of M and N that the probability Pr(n2|mi) in the first experiment is equal to the 
probability Pr(mi| ^2) in the second experiment. Therefore, a relation is indicated between the pairs representing the 
sequences A and A. 

For our purpose, it is sufficient to assume that the pair ^ that represents sequence A is determined by the pair a 
that represents sequence A in the limit as t2 ^ h, so that 

^ = R(a), (27) 

where the reciprocity operator R is invertible (since R(^) = a), and is assumed continuous. We shall assume that 
the above relation also holds more generally for sequences of arbitrary length. 

Now, consider the sequences A — [mi, 1712,1713] and B = [mi, 7172,1713], with 77x2 m'2, obtained from some experi- 
mental set-up, and the sequence, C, that combines these in parallel, namely 

C = A\/ B = [mi,{m2,m2),m3], (28) 

and take the limit as the times, ti, t2 and of the respective measurements coincide. The pair that represents C can 
be computed in two distinct ways, either as the pair R(c), or as the pair R(a) + R(b) that represents A V B . These 
two expressions must agree. Therefore, for any a and b, 

R(a + b) =R(a) + R(b), (29) 

which implies linearity of R, 

= ( ul f) ( ■ (30) 



-R3 R4 J \Q^2 

Similarly, by considering two sequences A and B that can be combined in series to yield C = A - B, and noting 
that ^ =B one obtains 

R(a0b) =R(b) R(a), (31) 
which, for any selected form of multiplication ©, constrains the reciprocity coefficients Ri, . . . ,Ri. 



3. Repeated Measurements 

Consider an experiment in which measurements M and N are performed at times ti and t2 respectively (see Fig. 5a). 
N allows only two atomic outcomes, 1 or 2. Sequences A = [m, 1] and B = [m, 2] have pairs a and b, respectively. 
Since either one outcome or the other occurs, P{A) + P{B) = 1, so that 

p(a)+p(b) = l. (32) 

Now consider an experiment where measurement M is performed at almost-coincident times ti and ^3, interleaved 
at intermediate time t2 by the trivial measurement N which has only one possible outcome (1,2) (see Fig. 5b). The 
sequence [m, (1, 2), m] can be written as 

C = [m, 1, m] V [m, 2, m] 
and, because the time-offsets are negligible, we also have that 



[to, 1, to] = [to, 1] • [1, to] = a • and [m, 2, to] = [to, 2] • [2, to] = B • ^ . 
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A B C 



(a) (b) 

FIG. 5: Experiment (a): For given m, the sequences A and B are mutually exclusive and exhaustive — the outcome of N must 
be 1 or 2, so that P{A) + P{B) = 1. Experiment (b): sequence C has associated probability P{C) = 1 in the limit as ti,t2 
and ts coincide. 

Therefore the pair that represents C is 

c= (a0R(a)) + (b0R(b)). (33) 

Now, the intermediate measurement, N, is trivial in that it only registers that a physical system is detected in the 
measuring device at time t2, but demonstrably does not affect the outcome probabilities of subsequent measurements 
performed upon the system. As measurement M is repeatable (see Sec. II), it follows that, in the limit as ti,t2 and ts 
coincide, 

p(a)+p(b) = l =^ P(c) = l. (34) 
We are now in a position to eliminate forms (C2) and (C3), leaving (CI) together with the specific form of p. 

Form (C2). Multiplication is via 7 = (1, 0, 0, 0; 0, 1, 1, 0), with p(x) = \xi\"e^='^/='K 

On substituting the linear form of Eq. (30) into Eq. (31), one finds that the only non-trivial reciprocity operator is 

R(.) -(;'). 

This is not invertible, which at once eliminates (C2). 

Form (C3). Multiplication is via 7 (1, 0, 0, 0; 0, 0, 0, 1), with p(x) = |a;i|"|x2|'^. 

On substituting the linear form of Eq. (30) into Eq. (31), one finds that the only non-trivial reciprocity operators 
are 

R(a)=(:j) «d R(a) =(:;). 

These are invertible, as required. 

Choosing R(a) = (a2, oi)^ makes c = {aia2 + 6162, cna2 + ^1^2)^ so that Eq. (34) reads 

|air|a2|'' + |&ir|b2|^ = 1 =^ |aia2 + bi62r+^ = l. 

The special case 6162 = — aia2 can satisfy the left condition while contradicting the right, thereby disproving this 
choice. 

The other choice is R(a) = (ai, 02)^, for which c = (af + 6f , + ^i)^' so that Eq. (34) reads 

|air|a2|'' + |6ir|62|^ = 1 =^ (a? + 6?)" {4 + bl)' = 1. 
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The special case ai = 62 = rt and 02 = 61 = r/t with r,t y^O reduces this to the identity 



valid for arbitrary t. This requires either a = 2 with /3 = 0, or a = with (5 = 2. But a = makes p(a) independent 
of ai, and /? makes p(a) independent of a^, whereas we require p{a) to depend on both arguments. Hence, this 
choice too is disproved, which ehminates (C3). 

Form (Gl). Multiplication is via 7 = (1, 0, 0, —1; 0, 1, 1, 0), with p(x) = (x\ -\- x^"'^'^ . On substituting the 
linear form of Eq. (30) into Eq. (31), one finds that the only non-trivial reciprocity operators are 



R(a) 



and 



R(a) 



ai 

-02 



These arc invcrtiblc, as required. 

Choosing R(a) — (ai, 02)^ makes c = (a^ — a| + 6f — 2aia2 + 26162)^) so that Eq. (34) reads 



(a? + ai)"/V(6? + 6r'^ = l 



{a\ - + 6? - 6^) V 4 [a^a^ + 6162)^ 



V2 



1. 



The special case 61 = 02, 62 = — oi can satisfy the left condition while contradicting the right, thereby disproving this 
choice. 

The other choice is R(a) = (ai, —02)^, for which c = (a^ + al + 6f + 621 0)^i so that Eq. (34) reads 



(a? + a^)"/%(6? + 6^)"'^ = l 



2\"/2 



(a? + a^+6? + 6i) =1. 



This requires a = 2, and this setting gives an admissible solution. Hence 

j)(x) = x\-\-x\. 

We are left with just this one solution. 



(35) 



C. Summary 



In order to combine sequences in parallel, we have the sum rule of Eq. (12), 




which we recognize as complex addition. In order to combine sequences in series, from Eq. (14) with 7 given by the 
surviving form (CI), we have 

which we recognize as complex multiplication. Hence the number pairs a, b, . . . behave as complex numbers, combining 
according to the rules of complex arithmetic. For the probability associated with a sequence, form (CI) gives Eq. (35), 
namely 

p{x) = xl + xl. 

These are Feynman's rules. 



IV. DISCUSSION 



In this paper, we have shown that the concept of complementarity, regarded by Bohr as one of the most fundamental 

lessons of quantum phenomena for our physical world-view [21, 26, 27], can be used to derive quantum theory. In 
particular, once complementarity is postulated in the minimalist form that there are two real degrees of freedom 
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associated with each sequence (or 'path'), but that one can only access one rcal-vahied function of these degrees of 
freedom in a given experiment, the complex arithmetic of the quantum formahsm emerges naturaUy. 

It is also interesting to consider which non-classical features the derivation does not use. It has been suggested [19] 
that the quantum formalism may owe at least a significant part of its structure to the fact that quantum theory permits 
non-locality and no-signaling to peacefully coexist, and a number of recent reconstructive approaches [12, 16, 28] rely 
upon postulates that concern the behavior of physically separated sub-systems. However, the derivation we present 
here takes place without making reference to more than one physical system, and thus demonstrates that features 
such as non-locality and no-signalling are not, in fact, essential to an understanding of the structure of the quantum 
formalism. That is, from the point of view of the present derivation, features such as non-locality and no-signalling 
are not fundamental, but secondary. 

The derivation also illuminates the nature of the relationship between the quantum formalism and the fundamental 
concepts of classical physics. The original formulations of quantum theory both make explicit use of the structure 
of classical physics Schroedinger's derivation [29] was based directly on de Broglic's wave-particle duality (itself 
based on the classical models of waves and particles), while Hcisenbcrg's derivation took the classical model of 
electromagnetic radiation from atoms as its point of departure [30]. Hence, both of these formulations presuppose 
the gamut of fundamental classical concepts such as space, time, matter, energy, and momentum. Therefore, the 
question naturally arises as to whether this is an historical accident and the quantum formalism is, in fact, prior to 
these concepts. 

In contrast to Schroedinger's and Heisenberg's formulations, the derivation presented here makes no explicit refer- 
ence to the classical concepts of space, energy, and momentum. The only aspect of the concept of time which has 
been assumed is the time-ordering of events; in particular, quantification of time (on the real number line) plays no 
role. Hence, the derivation shows that the core of the quantum formalism is a self-contained theoretical structure 
that makes minimal use of the fundamental concepts of classical physics. Of particular relevance to the programme of 
quantum gravity, the derivation clearly suggests that the Feynman rules are logically prior to the structure of space 
or, more generally, of spacetime. 

Since the formulation of quantum theory, numerous proposals have been made on how the formalism could be 
modified in various ways, such as allowing non-linear continuous transformations [31, 32] and modifying the formalism 
to use quaternions [33] orp-adic numbers [34]. Although these proposed modifications may be; mathematical plausible, 
they are rather ad hoc from a physical point of view since they are not clearly driven by physical facts or principles, 
and furthermore are difficult to subject to experimental tests. The derivation of Feynman's rules given here provides 
a natural framework within which such proposals can be systematically studied. 

Feynman's rules do not exhaust the content of the standard quantum formalism. For example, once translated into 
the von Neumann-Dirac state picture (see, for example [11, 20]), the Feynman rules imply that state evolution is linear, 
but do not imply that it is unitary. Therefore, one must appeal to additional arguments (such as Wigncr's theorem) 
to establish unitarity. It has been shown elsewhere that, given unitarity, the remaining standard structure — the 
tensor product rule, the representation of reproducible measurements by means of Hermitian operators, the general 
form of the temporal evolution operator, and the explicit forms of commonly used measurement operators — can all 
be systematically reconstructed [11, 13]. Hence, the present derivation provides a sound basis for the reconstruction 
of the entirety of the standard von Neumann-Dirac quantum formalism for finite-dimensional quantum systems. 
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Appendix A: Solution of ©-Associativity Equation 

Substituting Eq. (14) into Eq. (S3), and equating the first and second components, respectively, yields the following 
two equations: 

ill + 7375)ai6iCi + (71 72 + 7475)ai^'iC2 + (71 72 + 7376)^1 &2C1 + (7! + 7476)^16202 

+ (7173 + 7377)a26iCi + (7273 + 7477)^26102 + (7174 + 7378)026201 + (7274 + 7478)^26202 

= (7i + 7275)«i6iOi + (7172 + 7276)ai6iC2 + (7173 + 7277)^16201 + (7i74 + 7278)016202 

+ (7173 + 7475)026101 + (7273 + 7476)026102 + (73 + 7477)026201 + (7374 + 7473)026202 

and 

(7175 + 7577)01 61 ci + (71 76 + 7578)oi6iC2 + (7275 + 7677)016201 + (7276 + 7673)016202 

+ (7375 + 77)026101 + (7376 + 7778)o26iC2 + (7475 + 7773)026201 + (7476 + 71)026202 
= (7175 + 7576)oi6iCi + (7275 + 75)016102 + (7375 + 7677)oi62Ci + (7475 + 7573)016202 

+ (7177 + 7578)a26ici + (7277 + 7673)o26iC2 + (7377 + 7773)026201 + (7477 + 78)026202. 

These equations must hold for any 01,02,61,62,01,02- Equating coefScients, we get sixteen equations which, upon 
factorization and removal of redundant equations, reduce to twelve equations, 



7276 = 7475 


(Al) 


7377 = 7475 


(A2) 


74(72 - 73) = 


(A3) 


74(76 - 77) = 


(A4) 


75(72 - 73) = 


(A5) 


75(76 -77) = 


(A6) 


72(71 - 77) = 73(71 - 76) 


(A7) 


74(71 - 77) = 73(73 - 73) 


(A8) 


77(71 - 77) = 75(73 - 73) 


(A9) 


77(72-73) =76(73-73) 


(AlO) 


75(72 - 73) = 76(71 - 76) 


(All) 


72(72 - 78) = 74(71 - 76)- 


(A12) 



To solve the above equations, we select the nature of 76 and 77, choosing from the cases 76 = 77 7^ 0, or 76 ^ 77, 
or 76 = 77 = 0. 



1. Case 76 = 77 5^ 0. 

In this case, the first two equations give 

7475 

72 = 73 = , 

76 

while the remainder reduce to 

74(71 - 76) = 72(72 - 73) 
75(71 - 76) = 75(72 - 78), 

which both read 

, 75(72-78) 
71 = 76 H • 

76 

Therefore, 



f , 75/7475 \ 7475 7475 \ 

7 = 76 H 1 73 , , ,74; 75, 76, 76, 78 

y 76 V 76 / 76 76 / 
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which we can write in the more symmetric form 

7 = (0 - Ve, (j)e, (j); Oe, e,e,ip + (pe), (A) 

with real constants 6, (j), tp, e. 

2. Case ye ^ 77- 

In this case, Eqs. (A4) and (A6) give 74 = 75 = 0, and the remaining equations are 



7276 = 


(AlO 


7377 = 


(A2') 


72(71 -77) = 73(71 - 76) 


(A7') 


73(73 - 78) = 


(A8') 


77(71 - 77) = 


(A9') 


77(72 - 78) = 76(73 - 78) 


(AlO') 


76(71 - 76) = 


(All') 


72(72 -78) = 0. 


(A12') 



If both 76 and 77 are non-zero, then Eqs. (A9') and (All') imply 75 = 77 (= 71), contrary to assumption. Hence 
exactly one of them must be zero. Suppose, then, that 75 = 0, with 77 ^ 0. Then 73 = 0, 71 = 77, and 72 = 78, 
giving 

7 = (71, 72, 0, 0; 0, 0, 71, 72). (B) 
Similarly, suppose that 77 = 0, with 75 ^ 0. Then 72 = 0, 71 = 75, and 73 = 73, giving 

7=(7i,0, 73, 0; 0,71, 0,73). (C) 

3. Case 76 = 77 = 0. 

Before considering this case, we consider the solution of Eqs. (A1)-(A12) with respect to the nature of 72 and 73, 
choosing from 72 = 73 7^ 0, or 72 7^ 73, or 72 = 73 = 0. The treatment mirrors that of 75 and 77 just given. The first 
choice repeats solution (A) and the second repeats solutions (B) and (C). 

All that remains is 72 = 73 = 0, which we only need analyze in the context of 76 = 77 = 0. The surviving equations 
reduce to 

7174 = 7475 = 7578 = 0, 

whose solutions 

7= (71, 0, 0, 0; 75, 0, 0, 0) 
7 = (71, 0, 0, 0; 0, 0, 0, 78) 
7 = (0, 0, 0, 74; 0, 0, 0, 78) 

are special or limiting cases of solution (A). 

Hence, the possible solutions for 7 are the commutative solution (A) and the two non-commutative solutions (B) 
and (C). 

Appendix B: Solutions of Probability Equation 

We solve the probability equation, Eq. (26), for each of the five standard forms of 7 with the aid of two of Cauchy's 
standard functional equations 

f{xy) = f{x)f{y) and f{x + y)^f{x)f{y). 
We quote as needed [35] their continuous solutions, respectively 

fix) = and fix) = 
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Form (CI): 7 = (1, 0, 0, -1; 0, 1, 1, 0). 

Explicitly, Eq. (26) reads 

p{aibi - 0262,0162 + 0261) =p(ai, 02)^(61, 62) (Bl) 

for arbitrary oi , 02 , 61 , 62 . Change variables by setting Oi = r cos ^, 02 = r sin 0, 61 = s cos ^, 62 = s sin 0, with r, s > 0, 
to obtain 

p {rs cos{6 + 4>), rs sin{0 + 0)) — p{rcosO, r sin9) p{s cos <p, s sin </>). (B2) 
In case r = s = 1, this takes the form 

/(e + (/.) = /(^)/(<^) 

with /(V') = P (cos V', sin^), which has the solution /(V') = e^'^. Since f{ip + 27r) = f{il>), /3 = 0, so that 

/(■^) = p(cos V', sini/') = 1. 
Using this in Eq. (B2) with s = 1 and ^ = 0, we obtain 

p (r cos 0, r sin 0) = p{r, 0) , (B3) 

which reduces Eq. (B2) to 

p{rs,0) = p{r,0)p{s,0). 

This has solution p{t,0) = t". Hence, from Eq. (B3), p(r cos r sin <^) = r". Re- writing the arguments of p yields 

p{xuX2)^{xl+xiy^\ (B4) 
This satisfies Eq. (Bl), so is the general solution. 

Form (C2): 7 = (1, 0, 0, 0; 0, 1, 1, 0). 

Explicitly, Eq. (26) reads 

p(ai6i, 0162 + 0261) =p(ai, 02)^(61,62) (B5) 

for arbitrary oi, 02, 61, 62. In case Oi = 61 = 1, this reduces to p{l, 02 + 62) = p{l, 02)^(1, 62), whose solution is 

p(l,X2) = e''^^ (B6) 

In case 02 = 62 = 0, Eq. (B5) reduces to p(ai6i, 0) = p(ai, 0)p(6i, 0), whose solution is 

p(xi,0) = |xir. (B7) 

In case oi = 62 = 1, 02 = — 1/61 with 61 0, Eq. (B5) reduces to p(6i, 0) = p(l, — l/6i)p(6i, 1). Using Eq. (B6) and 
Eq. (B7), this gives 

p(6i,l) = |6ire^/*^ (B8) 
In case oi = 62 = 1, Eq. (B5) reduces to p{bi, 1 + 0261) = p{l, 02)^(61, 1). Using Eq. (B6) and Eq. (B8), this gives 

K61, 1 + 0261) = |6i|«e^(i+»=''i)/''i 
from which the solution can be read off as 

p(a;i,X2) = IxiPe^^^/^i. (B9) 
This satisfies Eq. (B5), so is the general solution. 
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Form (C3): 7 = (1, 0, 0, 0; 0, 0, 0, 1). 

Explicitly, Eq. (26) reads 

p(ai6i, 0262) = p(ai, 02) p{bi,b2) (BIO) 
for arbitrary ai, 02, 61, 62- In case 02 = &2 = 1, this reduces to 

p{aibi,l) = p{ai,l)p{bi,l), 

whose solution is 

p{xi, 1) = |a;i|". 

Similarly, by considering case ai = 61 = 1, we obtain p{l,X2) = \x2\^- Using these special solutions in Eq. (BIO) 
with (ai,a2) = (a;i,l) and (61,62) = (1,2:2) yields 

p{xuX2) = \xir\X2f. (Bll) 

This satisfies Eq. (BIO), so is the general solution. 

Form (Nl): 7 = (1, 0, 0, 0; 0, 1, 0, 0). 

Explicitly, Eq. (26) reads 

p{aibi,aib2) = p(ai, 02)^(61, 62) (B12) 

for arbitrary ai , a2 , 61 , 62 • The left side is independent of 02 , so p cannot depend on its second argument (02 on the 
right). Hence, p{xi,X2) — f{xi). Eq. (B12) thus reduces to /(ai6i) = /(ai)/(6i), whose solution is f{xi) = |a;i|". 
Hence, the solution of Eq. (B12) is 

p{xuX2) = \xi\''. (B13) 
This satisfies Eq. (B12), so is the general solution. 

Form (N2): 7 = (1, 0, 0, 0; 0, 0, 1, 0). 

Explicitly, Eq. (26) reads 

p{aibi,a2bi) = p(ai, 02)^(61, 62) (B14) 

for arbitrary ai, 02, 61, 62- Arguing as above, the left side is independent of b2, so p cannot depend on its second 
argument (62 on the right). Hence the solution of the above equation is also 

p{xi,X2) = \xi\". 

This satisfies Eq. (B14), so is the general solution. 
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